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^ ' Abstract 

^ : 

. The infrared behavior of the gluon and ghost propagators is studied in 517(2) Euclidean Yang-Mills theory 

||0 ' in the maximal Abelian gauge within the Gribov-Zwanziger framework. The nonperturbative effects associated 

with the Gribov copies and with the dimension two condensates are simultaneously encoded into a local and 
renormalizable Lagrangian. The resulting behavior turns out to be in good agreement with the lattice data. 
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1 Introduction 



$H ' The study of the infrared behavior of the gluon and ghost propagators has been the object of intensive investi- 
gations in recent years. Albeit not gauge invariant, these correlation functions enable us to probe the reliability 
of the various approaches which give rise to our current understanding of the behavior of Yang-Mills theories in 
the infrared, a task which is far from being achieved. This is due to the fact that propagators are the simplest 
Green's functions allowing us to evaluate in analytic form certain nonperturbative effects expected to be relevant in 
the infrared. Moreover, the lattice community has been able to develop accurate algorithms for a nonperturbative 
numerical study of the gluon and ghost propagators, which can be now analyzed on huge lattices, allowing therefore 
for a comparison between analytical and numerical results. Evidently, a qualitative agreement might be very en- 
couraging in pursuing further investigations of our theoretical frameworks. We also underline that this possibility is 
not restricted to a particular gauge. Nowadays, the gluon and ghost propagators can be studied from both theoret- 
ical and numerical viewpoints in several gauges as, for example, the Landau, Coulomb and maximal Abelian gauge. 



In this paper we focus on the gluon and ghost propagators in the maximal Abelian gauge [H HI 13] , which we 
are investigating since several years [H O [6l [71 [3 [9l [TOl [Til [121 O [H] ■ In our previous works we have provided 
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analytical evidence of nonperturbative effects which should be taken into account when facing the various features 
displayed by this gauge, such as the dual superconductivity picture for color confinement |15| and the Abelian 
dominance hypothesis [121 HZl [T51 (THj . It turns out that these nonperturbative effects can be accounted for by a 
set of dimension two operators which can be consistently introduced in the Yang-Mills action. 

As other Lorentz covariant gauges, the maximal Abelian gauge is plagued by the existence of Gribov copies pp] , 
requiring that the domain of integration in the Feynman path integral has to be suitably restricted to the so called 
Gribov region ^ 2T] . As discussed in [101 113| , this restriction can be implemented by introducing a dimension two 
nonlocal operator, known as the horizon function, namely 

5hor = d^xs-'A^ {M-'y'e-'^A^ , (1) 

where 7 is the Gribov parametei0 and (^A4~^^^^ is the inverse of Faddeev-Popov operator 

A yfab T-\ac T-\cb 2 ^ac J)d ac Ad /oA 

M = ~g e e A^A^ , (2) 

with Af^ and being the diagonal and off-diagonal components of the gauge field, respectively, i.e. A^ — 
and a = 1,2. Expression ([T]) generalizes to the maximal Abelian gauge the horizon function already obtained by 
Zwanziger [351 H3] in the Landau gauge. The nonlocal operator ([T]) can be localized by means of the introduction 
of a set of auxiliary fields (^", 0°, , ), so that the resulting action enjoys renormalizability [TUl [O] . 

In addition of the horizon function, other dimension two operators have been investigated. Our results have 
given support to the fact that the preferred vacuum state is that in which those operators condense, i.e. they 
develop a nonvanishing vacuum expectation value, lowering the vacuum energy of the theory. The first dimension 
two operator which has been studied is the gluon operator A'^A'^. This operator turns out to be multiplicatively 
renormalizable [24l [7] and its condensation, i.e. (^J^^JJ) ^ 0, gives rise to a dynamical mass generation for off- 
diagonal gluons [8], in agreement with the Abelian dominance hypothesis. As second example of dimension two 
operator let us quote the ghost operator e^^'c^c'', where c", denote the off-diagonal Faddeev-Popov ghosts. This 
operator is responsible for the spontaneous breaking of the global SL{2,M.) symmetry present in the ghost sector 
of the maximal Abelian gauge. It has been investigated by several authors [5S1 [S] |^ , see for instance ref.[T3] 
for a recent analysis of its renormalizability as well as of its condensation. The third dimension two operator which 
we shall consider is given by (0^0^ — i^fijjf — c'^c'^). It generalizes to the maximal Abelian gauge the operator 
introduced recently in the case of the Landau gauge ^28j. It reflects the nontrivial dynamics developed by the 
interacting auxiliary fields {4>f , (j>^ , cuf , ujf) needed to localize the horizon term ([1]). 

However, so far, these dimension two operators have not yet been analyzed simultaneously, a necessary step in 
order to get a more precise idea of their relevance in the infrared. This was due to the nontrivial task of explicitly 
constructing them, as in the case of the horizon function, eq.(IT]), as well as to the need of establishing their renor- 
malizability properties. The aim of this paper is that of filling this gap, by presenting a detailed analysis of the 
gluon and ghost propagators when all these dimension two operators are present in the starting action. In a sense, 
the present work can be seen as a kind of summary of our efforts towards a better understanding of the infrared 
behavior of the gluon and ghost propagators in the maximal Abelian gauge. 

The output of our results can be summarized as follows: 

• when all dimension two operators are simultaneously taken into account, the resulting local action remains 
renormalizable. This nontrivial feature is due to the large set of Ward identities which can be established 
when all operators are present. 

• the resulting behavior of the gluon and ghost propagators turns out to be in remarkable agreement with 
the available lattice data [551 133 131] ■ It is worth underlining that all dimension two operators affect the 
propagators. In other words, such a behavior of the gluon and ghost propagators can be obtained only when 
the dimension two operators are simultaneously taken into account. 

^We remind that the Gribov parameter 7 is not a free parameter. It is determined by the gap equation = 0, where F is the IPI 
quantum effective action evaluated by taking as starting point the Yang-Mills action with the addition of the horizon term | 22II23| . As 
such, the parameter 7 can be expressed in terms of the gauge coupling constant and of the invariant scale Aqcd ■ 
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The work is organized as follows. In order to provide a more easy reading of the paper, in Sect. 2 we give a general 
overview of our results about the gluon and ghost propagators, providing a comparison with the recent lattice data. 
In Sect. 3 we present a detailed discussion of the inclusion in the starting action of the aforementioned dimension 
two operators. In Sect. 4 we derive the set of Ward identities fulfilled by the complete action. In Sect. 5 we address 
the issue of the renormalizability of the model. Sect. 6 collects our conclusion. 



2 Summary of the results 



2.1 Notation 



Let us start by briefly reminding the standard notation in the case of the maximal Abelian gauge. The gauge field 
Afj. is decomposed as 

A^, - A^T^ = A^T^ + A^T^ . (3) 

The generator stands for the diagonal generator of the U{1) Cartan subalgebra of SU{2), while the index 
a — 1,2 labels the remaining off-diagonal generators {r"}. 

Accordingly, the field strength decomposes as 



f;, = DfAl-Dl^Al, 
Fl, = F^,=d^A,-d,A^+ge'''AlAl 

^ab ^3ab 



where we have introduced the covariant derivative with respect to the diagonal components of the gauge 
field, namely 

d;^" = <5"^a^ - ge^'A^ . (5) 



2.2 The tree level gluon and ghost propagators 



We collect here our results for the gluon and ghost propagators. 



• The off-diagonal gluon propagator: 

the transverse off-diagonal gluon propagator turns out to be of the Yukawa type 

m~k)Alik)) - (S,. - ^) S^" , (6) 

where m is the dynamical mass originating from the condensation of the gluon operator [S] 

= iA^^^ . (7) 

This behavior has been reported in lattice simulations [5U1 [ST] . It supports the Abelian dominance hy- 
pothesis, according to which the off-diagonal gluons should acquire a sufficiently large dynamical mass which 
decouple them at low energies. 



• The diagonal gluon propagator: 

for the diagonal gluon propagator we have obtained an infrared suppressed propagator of the Gribov-Stingl 
type, namely 
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where 7 is the Gribov parameter and is a mass parameter related to the condensation of the operator [32j 

Ojf = {4>^<t>t-u^t^t-c''c'')- (9) 

We observe that expression ([5]) does not vanish at the origin, in full agreement with the recent numerical 
data [31j . It gives rise to a positivity violating propagator in configuration space, a feature usually interpreted 
as evidence for gluon confinement. 

Moreover, it is worth to point out that the diagonal gluon propagator, fH]), can be naturally rewritten in 
terms of a power-law dynamical running mass of the type 

{A,{-k)A,{k)) = , (5,, ^\ , (10) 



where 

M\k) = -L^ . (11) 

Expression (jlOp is in accordance with the definition firstly envisaged in references [33l [34] and subsequently 
found in the operator product expansion (OPE) approach by 35 , and later in Schwinger-Dyson equations 
by [361. " 



• The symmetric ofF-diagonal ghost propagator: 

for the symmetric off-diagonal ghost propagator we have found 

i^i-k)c\k))„ ^ , (12) 

where z; is a mass parameter related to the condensation of the ghost operator |14j 

Oghost = ge'^'c'^c" . (13) 



Notice that expression (jl2p is suppressed in the infrared and attains a nonvanishing finite value at fc = 0. 
Again, this behavior agrees with that reported in [31]. 



• The antisymmetric ofF-diagonal ghost propagator: 

finally, for the antisymmetric off-diagonal ghost propagator we have 

(c"(-fc)c^(fc))antisymm = ^.4 + 2^i^k^ + (/i^ + y'^) ' ^^^^ 

As expected, this behavior is a consequence of the ghost condensate [14], {e°'^c°'c'') ~ v"^, being in agreement 
with [H] . 

In summary, the behavior shown above for the gluon and ghost propagators turns out to be in remarkable agreement 
with the most recent lattice data, as reported in [31] . This can be taken as a useful indication of the fact that 
the aforementioned dimension two operators play a relevant role in the infrared. Let also underline that all mass 
parameters, (m, 7, /z, i;), entering the gluon and ghost propagators are not free parameters, being determined in a 
dynamical way as solutions of gap equations, obtained by minimizing the vacuum energy, see for instance refs. [14l [H] 
for an estimate of the values of m and v at one-loop order. As such, all parameters will get proportional to the 
unique scale of the theory, i.e. (m,7,/x, u) oc A^qcd- 
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3 Identification of the complete classical action 
3.1 The Yang-Mills action and the gauge fixing term 

In order to obtain the complete classical action, let us start by specifying the gauge fixing term, namely 

So — Sym + 'S'mag : (15) 
where Sym is the Yang-Mills action in Euclidean spacetime 

SYM = \j A + F,,,F^,) , (16) 

with F^^,F^,y and £>^^ given in eqs.(l4|),([5]). The term S'mag in expression (fTSll stands for the gauge fixing term of 
the maximal Abelian gauge, being given by 

Smag = J A [ib'^DfA''^ - c'^M'^'^c'' + ge''''c''{Dl'A''^)c + ibd^A^ + cd^{d^c + ge^'^A^c'') ] , (17) 

where {b°',b) are the off-diagonal and diagonal Lagrange multipliers enforcing the gauge conditions, given by 
= and d^Af^ = 0. The fields {c"" , c"" , c, c) are the off-diagonal and diagonal Faddeev-Popov ghosts, 
respectively, and M.'^^ denotes the Faddeev-Popov operator 

^ab ^ -D'^^Df ~ g^e'"'e'"^Al^A'l . (18) 

The action is left invariant by the nilpotcnt BRST transformation 

= , (19) 



s^^ = -{Dfc>'+ge''''Alc), sA^ = -(<9^c + .ge'^^A^c^) 

sc" = ge'^^c^c, sc = f£°Vc^ 

sc" — ih'^ , sc — ib , 

s6° = 0, s6 = 0. 

Notice that the gauge fixing term p?]) can be written as a pure BRST variation 



(20) 



Smag = s J d^x (c'^DfA^ + cd^A^) . (21) 



3.2 Introduction of the horizon function, localization, and softly broken BRST in- 
variance 

As already mentioned, the maximal Abelian gauge is affected by the existence of Gribov copies, which have to be 
taken into account in order to properly quantize the theory. To deal with this problem it is necessary to restrict the 
domain of integration in the Feynman path integrals to the so-called Gribov region 17. In the case of the maximal 
Abelian gauge, this region is defined as (TUl US] 

n = {{A^, A^), DfAl = 0, d^A^ = 0, A^»^ = -Dl'Df - ^^^'^c^^'^A^Af, > } . (22) 

The restriction of the domain of integration is achieved through the introduction of the horizon function Shor, 
eq. ((!]). Therefore, for the partition function we write [TUinSl 

Z= /[rfA][d&][dc][dc]e~''^™+^"*<^+^'"') . (23) 
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The nonlocal term Shor can be localized by means of a pair of complex vector bosonic fields, {(ji'jf', <j)'^^) according 
to 

e-^'"' = j [d^] (det 7W)%xp I - y [ ^fM^^i^f + I'^ge"^ {<pf - ^f) ] | , (24) 

where the determinant, (detA^)^, takes into account the Jacobian arising from the integration over the fields 
(t>'^). This term can also be localized by means of a pair of complex vector anticommuting fields (lj^'', a)"''), 
namely 

Therefore, the horizon function gives place to a local term S'Locab namely 
e"^'"'- = y"[d(^][d(/)][dw][dcj]e-^^-'" , 

5Locai = j d\[ r^'M'^'^f - cS^^A^''^^ + j'gs^" {^f - r') ] . (26) 
Following [221 [23], we introduce the BRST transformations of the localizing fields {<i>'^ ,<i>'^) and (w°'',cD°'') as 

scj.f = , s< = 0, 

s< = <^ s^^ = 0. (^^^ 

It should be noted, however, that expression iSLocai does not exhibit BRST invariance, which turns out to be broken 
by soft terms proportional to the Gribov parameter 7. In fact 

/ j4_ 2 ab ( lab 7ab\ a 2 / i4 ^ [ ^cib ab a ^ab ( lab 7ab\ /a . ^mn Am n\1 /oo^ 

s d x"f ge [cl}^ - c/}^ ) ^ j g d x [e uj^ A,, - e [(pf, - (pf, ) [d^c + ge A^c )\ . (28) 



Nevertheless, as in the case of the Landau gauge [221 [231 [37] > the soft breaking ([28)) does not spoil the renormaliz- 
ability of the theory [TOl [13] ■ This remarkable feature relies on the possibility of extending to the maximal Abelian 
gauge the same procedure outlined by Zwanziger in the case of the Landau gauge [221 123j , amounting to embed 
'S'Locai into a more general action, S'lq^j^j, which enjoys exact BRST invariance, namely 

•S'Local — ^ "^Lo^al ; ■S'^Lj^cal — • (29) 

Furthermore, as it will be shown below, the term S'Locai can be easily recovered from S'l^^^j. The manifest BRST 
invariant action S'Locai is found to be [lOl [T3] 

d^x [r^M^^^^;^ - >(-c.f + ojfT-^<p';^ + M^iDi^cj^f + ^v^t^r^? 

+N;1[DI^u:!' + ge''^{d^c + ge''^Alc^)^^^] + M;1[DI^^^^ , (30) 

where 

j,ab ^ 2.g£'"=(9^c + ge'^''A'lc'')Df + ge''''d^{df,c + ge'^'^Af^c'^) - g\e'"'e'"^ + e'"^e'''')A'l{Dl^c^ + <?e""^;^c) , (31) 
and the external sources (M^^,M^^), {N'^l^N'^D transform as 

In order to reobtain S'Locai by the BRST invariant action iSl'^j^^i &rst take the physical limits of the external 
sources {M'^l,M'^l), [N'^l^N'^D, which are defined by [TUlin] 



^jab 



Arab 



phys 

phys ^'^ 



(33) 

= 0, 

phys 
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and then perform a shift in the variable o;"^ as [TUl [13] 

< - < + {M-'Y' [T^Uf + l^ge'\d^c + ge'^Alc')] , (34) 

so that 



Cfinv 
•^Local 



— •S'Local ■ (35) 

phys 



Thus, we consider the following action 

Si = So + S'Lo^'cal I (36) 

which enjoys the property of being BRST invariant, 

sSi=0. (37) 
3.3 Inclusion of the quartic ghost term 



Albeit BRST invariant, the action 5*1 is not yet the most general classical action to start with. The nonlinearity 

^ab A b 



of the gauge condition, D'^^A'' = 0, requires the introduction of a quartic term in the Faddeev-Popov ghost fields 



^ c^c^'^c'' , (38) 

which is in fact needed for renormalization purposes. In our case, due to the presence of the localizing fields 
(0^'', 0^'', a;|^^, w^''), the quartic ghost term is introduced in a BRST invariant way through the following invariant 
action Sa 

^ I ,A r_^.,« rib-a-b„ , 2-abiab f Jed icd ,-,cd, rd\ it lac —b b ' 



a 



-2z52,^-0-5''c'' + 2g3cD-0-e^''c''c^c] , (39) 

where a is a gauge parameter, which has to be set to zero after the renormalization procedure. In fact, introducing 
the action 5*2 as 

^2 - 5i + 5„ , (40) 
it follows that the equation of motion of the off-diagonal Lagrange multiplier b°- gets modified according to 



^ = iDfAl + a (5" + ige'^^'c^'c - ig'uJ^^<\}'^e) . (41) 



Therefore, one can see that the gauge condition of the maximal Abelian gauge, D'^^^'A''^ = 0, is attained in the limit 
a ^ 0, which has to be taken after the removal of the ultraviolet divergences. We also remark that the whole 
term Sa vanishes in the limit a — * 0, allowing us to integrate out the locafizing fields ((/)J^'', ^Jj'', wjj'', u)"''), and thus 
recovering the horizon function ([T]). 

3.4 The global U{8) symmetry 

In addition to the BRST invariance the action 5*2 displays a global U{8) [TUlIISj symmetry expressed by 

Q;IS2^0, (42) 

with 

Qf^ = f d'^xi (j)':"' -J-r - 4>t' — + cj"" -J-r - u)f -J— + — ^ - Mf, —L- + NZ — ^ - iV"^ ^ > 



J \^ 6<I)'=J' "^"^ %° Slu^'' "^'SM^'l, ""SM^l "^SN^l ^N^^ 



(43) 
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iV 


N 


dim 


1 


2 


2 





1 


1 


1 


1 2 


2 


2 


2 


gh. number 








-1 


1 








1 


-1 





1 


-1 


Qg-charge 














1 


-1 


1 


-1 1 


-1 


1 


-1 



Table 1: Quantum numbers of the fields and sources 



The presence of the global invariance i7(8) means that one can make use [T0l[T3] of the composite index i = {a, fi), 
i = 1, . . . , 8. Therefore, setting 

= , (44) 

and 

v-'"/ii/' -""/ii/i ^'/iiy' -"/ji// v^"/ju ^'^^/iii ^'/iu -"/ji^ I v^^y 

we can write S'2 as 

S, = 5ym + Smag + J d'x {^^X^'V^ - QIM'^'l^ + Co^T^"^': + M^^lJf + N^.^Dfco'l 

-2zg2<i,-0^-6fc'' + 2g3cSfC£^^c''c^c] . (46) 
For the symmetry generator we have 

C«^/A(*?A-*;^+.rA-.J^ + M;.jA__fl;;,_i^ . (47) 

By means of the trace of the operator Q,j, i.e., Qa = Qg, the i-valued fields turn out to possess an additional 
quantum number, displayed in the Table [l] together with the mass dimension and the ghost number. 



3.5 Introduction of external sources 



In order to establish the set of Ward identities, we have first to properly define the nonlinear transformations of 
the fields, as given in (PO)) . To this purpose, we notice that the BRST transformation of the gauge field can be 
written as the sum of two composite operators, namely 



where, 



/■r\ T^ab b 

Oi = -D„ c , 



O2 = -gs'^'Alc . 

Thanks to the fact that the BRST operator is nilpotent, i.e. = 0, it follows that 

sOi ^ ~s02 . 

These two operators can be defined by means of a suitable set of external sources, (ilj^, A'^, ^°), as 



-■ext 



To guarantee the BRST invariance of ^g^t we require that 



(48) 
(49) 

(50) 
(51) 
(52) 
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The nonlinear BRST transformations of the fields Af^jC"", c can be accounted for by the external sources ilf^,L°-, L, 
according to 



5, 



(2) 



s I d X 



where we require that 



Moreover, adding 5*^^^ to 5*2 we obtain an action that is left invariant by the following transformations: 



sf7. , sL° = , sL = . 



(53) 
(54) 



The 5i symmetry: 


The 5i symmetry: 




he" = wf, 















(55) 



As transformations ((55|) contain composite field operators, i.e., g£°-^4)\c and ge°'^u}^c, we define them by means of 
additional external sources {Y°-, Xf) and (Xf , Y"-), giving rise to two sets of BRST doublets 



so that 



^(3) 



sXf = 0, 

sy;" = 0, 



(56) 



ge Y^ (p^c-X, 



a ( ^ab b , 9 ^ab ^cd ib c c 



Ibc - ^e''''e'"^Lo'lc''c'^ 



+ge"^^Xtu>\c-Yt(ge"^\ 
Therefore, for the most general invariant external source term which can be added to 5*2, we obtain 

where the last term, which can be written as an exact BRST variation 

X I d'x (Af«M« +7v;,iv^) - -xs I d\N;,M;, , 



(57) 



(58) 



(59) 



is allowed by power counting and has to be added for renormalization purposes. Also, the parameter x stands for 
a free coefficient. 



3.6 Introduction of dimension two operators 

The last step towards the construction of the complete starting classical action is the introduction of the three 
dimension two operators Oa'^iOjj and Oghost, eqs.([71).([5|l. p^ . Let us start by considering the gluon operator 
0^2(2;) = i ^J^(x)AJ^(a;). Introducing the BRST doublet of sources (A, J) as 

sA = J , s J = , (60) 

it turns out that Oa^ can be introduced in a BRST invariant way, namely 

Sj^s f d^xX {Oa2 + ^CJ)^ f d\ {JOa2 + ^CJ^ + \ AlDfc"") , (61) 
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where C is a constant parameter needed to account for the uhraviolet divergences of the vacuum correlation function 

{{Al{x)Al{x)){Al{y)Al{v))). 

The operators O/^jOghost can be introduced in a similar way. More specifically, defining the BRST doublet 
of sources (r, a) as 

ST = a , sa — , (62) 
the invariant term S^y describing the coupling of Ojj is given by 

Sa = s j ct^xT {Off + CT + p J) = j dt^x {aOff + ^Ka^ + paJ - T sOjf) , (63) 

where k and p are constant parameters, needed for renormalization purposes. Notice in fact that expression 
contains the mixing term aJ. This term, allowed by power counting, accounts for the ultraviolet divergences of 
the mixed vacuum correlation function ((A°(a:)A^(x))(0^(?/)(/)^(y) — u}\{y)uj\{y) — c''{y)c^{y))). 

Finally, the introduction of a third doublet of sources (ry, 9) 

sri = e , se^O , (64) 
allows us to introduce the ghost operator Oghost(a;) = ge'^^c°'{x)c''{x), namely 

Sg^s j d^xr^ (Oghost + \I30)= J d^x (SOghost + i/3 0^ _ ^ sOgho^t) , (65) 

where /3 is a constant parameter needed for the divergences of the correlation function ( {e^^c" {x)c^ {x) ) (e™"c™ (?/)c" {y))). 
Notice, however, that the ghost operator Oghost breaks the symmetries (I55p . Therefore, to maintain the symmetry 
content of the theory is necessary to introduce two more BRST doubles of external sources. 



srii = ~9i , sOi = , 

s9i = fji, sfji = 0, 

and define an extra term given by 



(66) 



5extra = S J d^X gs'"' {O^cj)^ c" - 7J,iU^ c") 

d'-x [ge^'' (fy.^^c'' + ?7,0^c^ + O^uj^c'' + O^Cj^c'') - g'^e.^.c^c + g^m^tC'c] . (67) 



3.7 The complete classical action 

We are now ready to write down the complete classical action E, given by 

^ — So + ^cxtra 5 

50 = 5*2 + 5*0x1 + Sj + Sa + Sg , 
S2 = Si + Sa , 

51 = So + S'lo^j^j , 

5*0 = Sym + 5'mag , (68) 

where, Sym, 5'mag, S'^^^^i, S^, 5cxt, Sj, S^, Sg, 5cxtra are given, respectively, by ((111), ((TT]),®,®, (El, (Ell), 
(|67p . Thus, the complete classical action is 

S = 5yM + 5mAG + ^Lo^cal + 5q + 5cxt + Sj + Sa + Sg + 5cxtra , (69) 
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or, explicitly, we have 



YM 



a \ Aah lb 



fit r^t 



^ I ^2—aia/lbib —b, ,b\ o^2-,a /a-fc„6 



+s j d^xX {\AIA1 + lQj)+s j d^x T [{4>t<i^t - Cjfu:'} - c'^c'^) + lno + pj] 
+s J d'^x [ge""^ {tjc'^c'' + 9,(j)'}c^ - ri.Lolc^) + 

Sym + J d\ [ib'^DfAl - c'^W'^c'' + ge''''c''{Dl^Al)c + ibd^A^, + cd^{d^c + ge^'A^c'') + .^^A^-^V- 

+M-^ADt4>\ + ge'-'id.c + gs'''A'^^c'')Q\ ] - ^^Dfc' - ge'^'K^^Alc + \ ge'^'iDl^c^c 

c-Xf 



ab cdj^b c d 
2 M 

,y_^ab^cd^b^c^d 

a , 



n^id^c + gs'^'^A'^^c'') + ge''''L''c''c + |e°^Lc'^c'' + ge''^"'^''' 



^ab va—b^ 



ge^'u^lc 



+2g^Cjt€ e^^c^c^c] + \J A^A^ + XA1Dfc^ + a (^0" - ^^t^t " c°c°) + r (i6°c° - ge'^^c'^c^c) 

+ ^ J2 + p + ^ g^ab l^g-a^b ^ -^^a^b _^ ^^^a^b ^ g^^a^b ^ Q^^d'^C^) - igs'^Sb'^C^ - g^rjC^C^C 

2^ j.a a I 2_ -a a , n2\ 

-g f^iVt c c + g r]iUJi c c + - j . 



(70) 



The expressions for the gluon and ghost propagators given in Sect. 2 are easily derived by considering the relevant 
quadratic terms of (|70p and by replacing (J, a, 6) by the more conventional mass parameters (m? , /i^, w^) originat- 
ing from the corresponding dimension two condensates, i.e. m? ~ t8j, ~ {{4>i4'i ~ ^f^f ~ c'^c'^)) 132, 



(e c c 



The action ([70]) constitutes our starting point in order to establish the renormalizability of the model. 



4 Ward identities 



It turns out that E fulfills the following set of Ward identities: 



• The Slavnov- Taylor identity: 



5(S) 



d% 



6Y. 



ST, \ (51] 6^ ST, ST ST 



ST 



SK-J SA- 

•' L 



/i2 



ST ST ,ST ST 

-n— T y TT^ + 7 1- 

(5i'^ Sc"^ SL Sc Sc°- Sc 

^ ^ ""'s^- ^ sx- 



-Xf 



ST 



ST ST ST ST _ ST 
SYt ' ^^+^57+%"^'^+'^^^ 



= 



(71) 



The four global W^^^'' -identities which mix the Faddev-Popov ghost fields with the auxiliary localizing fields: 



w\^\t)^Q, (iV = 1,2,3,4) 



(72) 
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n 


K 


? 


L 


Y 


Y 


X 


X 


A 


r 


J 


O" 


V 


e 


Vi 


?7» Si 


0^ 


dim 


3 


3 


3 


4 


3 


3 


3 


3 


2 


2 


2 


2 


2 


2 


3 


3 3 


3 


gh. number 




-1 


-2 


-2 


-1 


-1 





-2 


-1 


-1 








-1 





-1 


-1 


-2 


Qg-charge 














1 


-1 


1 


-1 




















1 


-1 1 


-1 



where 



Table 2: Quantum numbers of the external sources 



6J: (5S 



Sc'^ 



(5E (5S 



(51] 



(5wf 

. (5E (51] 



/ (5S 



(56*, (5L 



<5S] 



(5E (5S] 



,51] -„ ,5S „ 

+^uii-^ - ^mTT^ + 0-, 2r?i — 

^',5^° ^*(5f^^ 6r,, SL 



5L^ 5Xf 
ST. 



6T ST 
-e— + 2r^,— 
orji oL 



ST 
SuJf 



' ^ 5L-^ Sij? 



" ,5L« S6? 



(73) 



The global U{8) invariance: 
with 



Q.,(S])=0 



(74) 



- ui- to- 

' Sujf J Sijf 

J * 



" SMZ 



'SN: 



HI 



' SY'' ^ 6Y^ 



' SX^ ' SXf 



m-c ?7jT^ 



69, 



(75) 



The trace of ([75]) defines a new charge displayed in the Tables [T] and [H This operator generalizes that of 
eq.p7|. 



• The exact rigid symmetries: 

m,,{T) EE 



5n(i'(s) 

9=l(2)(I]) 



SM'' 



ST 



9-^ 
'Sv, 

ST 

S^ 



0, 



'■SN 



SXf 



ST 



Sujf ^''SMf, ^'SNf, ' SY.r 



- ST ST\ ^ 



SY/" 



Yf" 



ST 
5X? 



(76) 



The diagonal gauge fixing condition: 



— - lOf^Af, 



(77) 
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The diagonal anti-ghost equation: 



Sc ' '"6n, 



(78) 



The SL{2,R) symmetry: 



I>(E)./a(c-§ 



<5E (5E 5Y. „,(5E , „ 
261— = 

SL" 5b'' 6L 



(79) 



The local U{1) invariance: 



with 



>V3(E) = -id^, 



6 A, 



ye Off 



(80) 
(81) 



where 



Off- 



{A;,b^, c", c^ 0^ 0^ c.?, 1^;:, i^;, e^, L^ , , m;„ m^"„ iv» , iv;,} . (82) 



• The BRST on-shell invariance of the general operator {0^2 + aOj^): 

^(E).y.^.(-+.--.c'^--2,-j.O. (83) 



5 Algebraic characterization of the most general counterterm 



We can face now the issue of the renormalizability of the starting action E. We shall employ the algebraic 
renormalization [38j and look for the most general invariant counterterm which can be freely added to all orders 
of perturbation theory. To that purpose we perturb the classical action E by adding an arbitrary integrated local 
polynomial E^x in the fields and external sources of dimension bounded by four, zero ghost number and zero 
Qs-charge. Requiring thus the perturbed action, E + eEcT, satisfies the same Ward identities as E to the first 
order in the perturbation parameter e, we get: 



5(EH 


h cEct) 


= + 0(6^), 




V cEct) 


= za^A^ + 0(e2), 




h cEct) 


= + O(e2), 


(EH 


h eEcx) 


= + O(e2), 


9^y(EH 


V eEcx) 


= + O(e2), 


^n(^)(SH 


h cEct) 


= + O(e2), 


Q,j(EH 


V eEcx) 


= + O(e2), 


I?(EH 


V eEct) 


= + O(e2), 


W3(EH 


V sEct) 




W(EH 


V cEct) 


- + O(e2). 
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This amounts to imposing the following conditions on Eqt 



5_ 

Tc 



,.,S(JV) 

KVj 2jCT = 

lA ScT — 



where iSs is the nilpotent linearized Slavnov- Taylor operator, 

= 



[N. 



1,2,3,4) 
1,2) , 



(85) 
(86) 



5T 5 



6T 



ST 



Sn- SK-) SA- 
S_TS_ 5TS_ 
Sc" SL'^ SLSc^ Sc SL ' 

A A 
V-'^^a 



ST 

\sn- 



SK- 



ST 6 



Sc"- Sc ' S(h^ * Sujf 



5T S 
SAf, SQf, 



""SM":. 



ST S 
SL" Sc"- 
S 



5N"^ 



SX" 



<5A 



St 



Srj Srji 



se,. 



while yy. 



(87) 



with N — 1, ... ,4, and 2?s are the linearized operators corresponding to the Ward identities (|73 



and (|79p . respectively, and they are given by 



S(l) 



S(2) 



S(3) 



d^a; ( 0,"-^ + + M",J- 

^^'Sc" S(j)f ^'sni 

L.c"—-N".— -X" 

Sivf '''SVLl ' SL" 

S ST S ST S 
+ 1- 



SL" 
5 



Sc" 
ST 

SY" ' 



S(4) 



ST S 
''SL"~Sdf} 

d^x 



c"—- M" — 

sujf ^'se 



Sc" Sc"SYf ' "SX^' 
S 



ST S .ST 


S 


SY^ 5b" "Sb" 


SYt 


ST S .ST 


S 


SX^ Sb" "Sb" 


SXf 


.ST S 1 


^ ST 


^'Sb" SXf ^ I 


yS^t 



ST 
SX" 



ST S ST S ST 5 

— i 



Sc" Sc" SX" 



SY" Sb" 



— — - + /V" — - AT" — + 0— - 2- — 

'SL"Sijf S4>t ^M^^ ^',5^!^ S-q, ^'^T. 



Sb" SY," 
5_ 
SL 



ST 
Suj" 



and 



Vy^ ^ d\ c" — 



ST S ■ 5 S 

^WSL"^ SL 



Sc" SL" Sb" 

For further use, let us write some useful commutation and anticommutation relations 



} 



S(3) 



{IHy, 5e} = Qij , 



S(4) 



9^(2) 



^_s_ ^ 2^ _S_ 

Sfji SL 



29, 



X" 



_S_ 

Sdi ' '"'SL 
5 
SL" 



— Y" 



S 
SL" 



(89) 



(90) 



From the second and the third constraints of ([85)l it follows that Sct is independent from the diagonal Lagrange 
multiplier 6, and that the diagonal antighost c enters only through the combination {Q^ + 9pc). Furthermore, 
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from general results on the cohomology of gauge theories [33, it turns out that the most general solution of the 
constraint iS^ScT = 0, i.e. the first of eqs. (|85p . can be written as 



ScT = ao 5ym + 5eA(-i) , (91) 
with A*^"^^ being an integrated local polynomial with ghost number —1, given by 

+a7 + d^,c)A,, + as c"i" + ag cL + ge''%\c + an X>," + ai2 ge''''Q\c 

+ai3 Yt^l + ai4 Ytr, + ais X^Q'^ + aie N% d^c^l + aiy N% ge'^'A^^I + a.gM;^ 9^cDf 
+ai9 m;, ge^'A^Oj', + aso ib^c^ + asi ge^'c^c'c + a22 ge^'c^A^Al + aaa 

+ a24 io'^cj.fojy^ + a25 cSr0'0"0? + «26 + 027 CS,''^?^^^" + ^28 ^f^^cD^'^^" 

+a29 cSf'/'j ^''/'j + 030 + a3icD,''(/.^"cS,^c.^^ + a,2 ^I'K^^c^ + 033 c^^^^c'^c" 

+a34 cS>.'c''c'^ + aas (^1¥IA^A^ + age cS^^^t^Jl + 037 ^1<\>\AIA\ + agg cSf^'C 

+a39 ge^'A^a^^^ + a4o ge''\dM<i>\ + ^41 X^V^Af/J, + i(a42 A + a43 t)AIAI 

+ i(a44 A + a45 t)A^A^ + (a46 A + a47 r)c"c" + (048 A + a49 + (ago A + 051 t)<I.>," 

+ (052 J + 053 cr)u)°;(j)°; + (a54 J + 055 (t)A + (age J + 057 (t)t + agg 77 9^ + agg 7/ 61 

+g£°'' (ago 11 c'^c^ + aei 77 C'/'^ + 062 ?7 + ^63 ^ w^^, + 064 e,(j)'lc^ + agg c'')] , (92) 

where the coefficients a„, n = 0, . . . ,65, are free dimensionless parameters. Notice also that in the derivation of 
expression (|92|l use has been made of the fact that the action E, and thus Sct, are left invariant by the following 
discrete symmetry 

y^~^y\ y^^-y"", ;ydiag ^ _^diag^ 3^ ^3^, (93) 

where 3^°, with a = 1,2, stands for the elements of the off-diagonal set ([5^ . while 3^<^'^s for the diagonal sector 

e [A^,h,c,c,n^,L,rJ,e,^l,,f,^,0^,9^] , (94) 
and 3^ the sources A, r, J, a. As one can easily recognize, this symmetry plays the role of the charge conjugation. 

After a quite lengthy calculation one finds that the most general expression for A^^^^ compatible with all constraints 
([55]) and the discrete symmetry ([M]) is 

A^-^) = j d%{{a, + a2){n'lA'l + ge'^%Alc) + [a, + a,)C,Dfc' + (ai - a4 + a^){N;,Dfct>\ - M^^Dfco^) 

-(ag + aii)ujtM''U\ - 06 c'^DfAl + L^c" + agLc + an x N;^^M;^, + a2o(^c"&'^ - ge''''c''c''c) 
-(-aas - a2o)g^tof<Pn4>';^' " + (-"Og - 2a2o)g^Lut^^c''c'' + i(o8 - og - aa43)A 

+ ia43 tA^^A"^ - (-aag + 2a2o - 0047) A(0,>," - uj^uj^ - c^c") - (ag + 047) t(0^9!'° - + 047 tc"c" 

-(ag - aQ)ge"-''{'q^c^ + 9i(f)fc'' - Ty^w^c^) + a54 JA + agg ctA + age Jr + a57 err + agg 776*1 . (95) 



Noticing that 

and renaming the coefficients as 



Jt ^ aX + S^iXr) , (96) 



01+02^01, a2 + a4^-a2, ag ^ 03, 

ag ^ a4, ag ^ as, a4i ^ ag, 020^— 07, 

1 /■ 

043 — * og, a47 ag, a54 aio 2^,7 
055 + 05g ^ an p, 057 ^ ai2 5K, ogg ^ ai3 (97) 
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we get 
A(-i) 



2V"3 ^"4 



-(as - a^)ge°-^{'qc''c'' + 6*^0^ c'' - r?,w"c'') + 5C -'^A + an pcrA + ai2 ^Kcrr + ais \f3i^e^ 



(98) 



so that for the counterterm ScT we get 

n 



^CT 



9 



d X \ (ao + 2ai) 

,2 



(A A A'^ A"- 4- A A A'^ A'^^ 



+ ^ (a^A,)(a^A, - a.A^J + (ao + Aa^)^-AlAlAlAl 
-aQi^a _ c«ge«''(a^A^)c'' - 2c"ge"^yl^5^c^ - g^c'^eA^.A 



\9 Aa Aa Ab Ah 



+i(ai - a4)6° 4 - (as + 04) 



+ (as + ai) 



C^'C - 0,?ge«''(a^A^)0,^ - 2^1ge'^'A,d^cty\ - g^^cf^fA^A, 



-(as + a4) [iold^u^ - u^gs''''idM^' " 2u,^gs''''A^d^u'^ - g^cofuj^A^A 
+ (2ai - as - a4)gh^'e'"'ic^c'' - 4>t<l^\ + Oj^^DA'^^A^ + (ai - 04 - a5).g£"''cVi?;^^^^ 
-(ai - as)(r!^ + d^c)ge-''Alc' + {a, - 2as - 04) We'^'e^^Co'^ A';,c''d^<p\ + 5'£"''£^'cD,r9(A^c'*)/ 



(2ai - as - 04 - a5).g3(,5°'=e'"^ + 



xu>tAlAl4>\c ~ (a2 + 2as)g2e-^£-''(7V^(/)^ + M^^CoDAIc^ + (a, + a, + a^^W^^D^ 



a^KlDfc" 



+a^ge''''KlAlc - {a, + + + a^)ge'''C^{Dl-c'')c + (as + 2a:if- e'^h^" C^Alc^ - a^ge-'L^c'c 



-{2a,-a,)^e''''Lc- 



'u;'i{d^c)cj,\ + 2g^0ot<t>1A^d^c 
M^^LoD - (ai + a2 + a^){M;,Dfcl>\ + N^^Dfu^ + M;,Df-^\ 



+X^oj\ - Yt4>\)c - (as + a4 + 05) \2ge''^u:'^{d^c)d^cp\ + ge''''u:'^(d^ c)4>\ + 2g^ujt^'i A^d^c 
+ (ai+a2 + a3 + a5)5£"'(a^c)(7V;,(/.^ ^-'^ --"^ '~ — , ^ ^/ ^ab .6 . a ^ab, 

+iV;,Z)f - as X(M« + N^.N^,) + ^ [a^fe^fo" + 2z(a7 - a^)ge'^%'^c'^c - g\a, - 2a^){^1^1 

-uIujI - -ec'')[^<i>] - ~ c^'c'') - 2^5^(07 _ 2a3M^1b''c'' + 2g\aj - 2as - 05)^^ (/.°e''"c^c^c 

+ (2ai + as - 04 - aafi)\j A^'^A^ + (ai - 04 - aa8)AAj;^i:>);^c^ + (03 + ag){ac°'c'' ~ irb^'c'') 

+ (as + as + a9)T5e''^c''c''c - (2as - a5)ge"'\eif'c^ + f],(l)'^c'' + n.^fc'' + e.uj^c'' + 9,cj^c'' - n^b^-c^) 

+2g2a3 (rycV'^c + ^.(/.^c'^c - ry,cD,V°c) + agiaA^A^; + a^TAlDfc' + a[a^ + ay + ag) 



-cSf cj," - c°c°) + A(i5'^c'^ - ge''^(fc^c) - (as + a9)(T(C0" - + aio^ + an p J(t 



hai2- cr 



/3 



(99) 



After the characterization of the most general local counterterm Eqt compatible with all constraints, eas.(|85p. we 
still have to check if it can be reabsorbed through a multiplicative redefinition of the fields, sources and parameters 
of the starting action E, according to 



E[*o, ^0, Jo, ??o, f^o, ^0, Ao, To, Jo, fio] = S[*, V, J, t?, f^, i^. A, r, J, ct] + e Ect + 0{e^) , 



(100) 



where 



-00 



Jo = Zj J, 



(101) 
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with 



J ^ {^;, L^ L, n^, , X", r/\ F,«, m;„ m^^, n;„ n;„ jj, e, m, h e.} 

= {5, a, X, C, /?} • 



(102) 



Moreover, by taking into account the mixing of the sources displaying the same quantum numbers, i.e, (fij^, -ft^^), 
(A, r) and (J, a), we shall set 



O " 



where the Z-matrices are given by 
^QK = I + e 



Z 



Ja 



Z\ 

ZtX Zt 



Zj Zja 
ZcfJ Z(j 



By direct inspection of Sct, the renormalization factors are found to be 



,1/2 



1/2 



_ -7 7-1/251/2 
2-1/2 ^ 2-1/2 

^ 



yl/2 _ 7I/2 

'5 — , 

7I/2 _ 7 

„ 7I/2 



7I/2 _ 7 7-1/2 7I/2 



z 



Zx = Z. 



N — ^g^c 
1/2 



ry _ /v — 1 7I/2 7 

— ^g ^M, 



7-2 7 7-1/2 
Zj„ ZjfiZjc. , 



7 7_ '7—1 '7I/2 7—1/2 

■^F — -^Y — -^c , 



Zn 



y^l 7I/2 7-1/2 

, 



7 7_ 7 — 2 7 7I/2 



Za 



-1/2 



7-171/27-1/2 

Z^ = Zmi 



z^ = z-^zli^z-'i\ 



Zn 



7I/2 



(103) 



(104) 



(105) 



with 



7I/2 

zl'' 
zl'^ 

Zm 

Za = 

Z^ = 
Zp = 

z^ = 

Z0 = 



= 1 + e 

= 1 - e 



ai 



2 

as + a4 

2 ' 
03 - 04 _ 
2 

ai + a2 + 



Zg = 1-e 



05) , 
as - ttA 



= 1 



T ' 

as - 04 



) 



1 + e(2ai + 2a2 + as - 04 - ae) , 
1 + e(ao + 2a4 + 07) , 

2ao — 2^1 + a^^as + 2a4 — 2a^a-j + 2aa8 — 2a^aQ + aio 



1 + e 

1 + e 



ao — (1 + a)a3 + 04 — aar — ^- 

1 — e 2a4 + 2— ag — 2a9 — 012 , 
L K J 

1 - e(ao - 2a3 + 2a4 + 2a5 - ais) , 



c 

ag) + an 



(106) 
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and 



= I + £ 





02 





-fli - ^(03 - a4) 



z 



Jo- 



-i(ao - 03 + 04) - aas 


a-s 


+a{a3 + aj + ag) 




i(ao - 03 + 804) - ag 


— oq + 03 — 04 — aas 




as ^ 



+a(a3 + fly + ag) 



(107) 



04 — Og 



This concludes the proof of the renormaUzability of the complete calssical starting action E. 



6 Conclusion 



In this paper the gluon and ghost propagators have been investigated by taking into account the effects of the Gribov 
copies as well as of dimension two operators. The output of our results is summarized in Sect. 2, where the ex- 
pressions for the tree level propagators can be found, being in good agreement with the most recent lattice data [31] . 

Certainly, much work is needed in order to reach a better understanding of the maximal Abelian gauge. Nev- 
ertheless, the results which we have obtained enable us to strengthen the fact that the agreement with the lattice 
data has been obtained only when the effects of the Gribov copies and of the dimension two operators have been 
simultaneously encoded in the starting Lagrangian, which enjoys the important property of being renormalizable. 
This point can be better clarified by the following considerations: 



• The quantization procedure and the issue of the Gribov copies 

The starting point to analyze Yang-Mills theories at the quantum level is by means of the Faddeev-Popov 
quantization formula, based on the introduction of a gauge fixing and of the corresponding ghost term. It is 
known that such a procedure is plagued by the existence of the Gribov copies. A full resolution of this issue, 
amounting to restrict the domain of integration in the Fcynman path integral to the fundamental modular 
region, is still unavailable. A partial solution to this problem consists of restricting the domain of integration 
to the Gribov region il, which is still affected by Gribov copies. Although this procedure does not eliminate 
all copies, it has the advantage of being effectively implementable. As we learn from the work of Zwanziger 
[22I [23^ in the Landau gauge, the restriction to the region is achieved through the introduction in the 
Yang-Mills action of a nonlocal operator, known as the horizon function. This nonlocal operator can be cast 
in local form by introducing a set of additional localizing fields. Remarkably, the resulting local action turns 
out to be renormalizable jUl [53] . This procedure has been successfully adapted to the maximal Abelian 
gauge [TOl [13] . A second point to be noticed is that the introduction of the horizon function in its local form 
is equivalent to the introduction of a specific dimension two operator. In fact, the Gribov-Zwanziger gap 
equation |2H 1221 123] determining the Gribov parameter 7, namely 

^=0, (108) 
07 

with r being the IPI effective action, is equivalent to require the existence of a nonvanishing dimension two 
condensate. In the case of the Landau gauge, this condensate is given by [2H H3] 

{f^^^A^ix){cb^^{x) - $^^{x))) ^ , (109) 

where (j>^^ , (f>^'~^ stand for the localizing fields and the indices A, B, C belong to the adjoint representation of 
SU{2). In the case of the maximal Abelian gauge the corresponding condensate is given by (e'^^A^(a;)(0°^(a;) — 
(l)'^{x)). The same feature holds in the Coulomb gauge, see [39] for a review. 

• Introduction of the dimension two operators 

As mentioned before, the inclusion of the horizon function is equivalent to the introduction of a dimension 
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two field operator in the localizing fields. Therefore, we can look for other dimension two operators which can 
be added to the theory, provided one is able to maintain renormalizability. From this point of view, the intro- 
duction of the three dimension two operators Oa^ = \ A^'^A"^, O// = {(t>\4>\ - - c'^c''), Oghost = (e^^c^c^) 
looks very natural. It is remarkable that these three operators can be simultaneously added to the horizon 
term in a way which preserves renormalizability. We also notice that all three operators considered here have 
their analogue in the Landau gauge, see [71 128j and refs. therein. In much the same way as the horizon 
function, these operators carry nonperturbative information, encoded in the corresponding condensates. 

The good agreement of our results with the lattice data can be taken as evidence of the fact these di- 
mension two operators play a relevant role in the infrared. For example, without the introduction of the two 
operators O//, Oghost, the infrared behavior of the off-diagonal ghost propagator would be deeply different 
from that of eg. (fT2|) . Instead, it would have displayed an enhanced behavior of the type as reported 

in our previous investigation JOj , where only the horizon function and the gluon condensate 0^2 were taken 
into account. The same occurs for the diagonal gluon propagator, eq.®. Without the introduction of Ojj 
it would be vanishing at fc = 0. 

We remark that the same features have been detected in the Landau gauge, where the most recent lat- 
tice data [40 l [4 H |42] point towards a finite and nonvanishing gluon propagator at fc = 0, while exhibiting a 
less enhanced ghost propagator. As discussed in ^28j these features can be accounted for by considering the 
effects of dimension two operators, which nicely fit within the Gribov-Zwanziger framework. 
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